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ABSTRACT 


In this study, the wake shape under symmetrical flight 
conditions and its effects on aerodynamic characteristics are 
examined. In addition, the effect of wake shape in sideslip and 
discrete vortices such as strake or forebody vortex on lateral 
characteristics is presented. The present numerical method for 
airplane configurations, which is based on discretization of the 
vortex sheet into vortex segments, verified the symmetrical and 
asymmetrical roll-up process of the trailing vortices. Also the 
effect of wing wake on tail planes is calculated. It is concluded 
that at high lift the assumption of flat wake for longitudinal and 
lateral-directional characteristics should be reexamined. 
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1 . INTRODUCTION 


In the past, considerable interest has been shown in 
computational methods of analyzing the roll-up characteristics of 
aircraft trailing vortices and their effect on tail surfaces. 
Although the general viscous flow problem is quite complicated and 
requires further research, it is possible to calculate the main 
effect of wake roll-up by a potential flow model. Considerable 
progress has been made in the latter. 

The occurrence of a vortex sheet which trails from a finite 
aspect-ratio wing is a consequence of the nonuniform lift 
distribution on the wing sections. The roll-up of this sheet into 
discrete vortices, as a result of a convective motion, was 
recognized by Lanchester in 1907 (Reference 1). 

Westwater (1935) modeled the continuous vortex sheet (assumed 
two-dimensional) trailing from a wing with an elliptical load 
distribution by discrete vortices, and computed the time development 
of the configuration by use of a complex potential. He assumed that 
the sheet extended to infinity upstream and downstream of the wing; 
and the effect of the wing itself was ignored. That is, no attempt 
was made to incorporate the effect of the solid body in the flow. 
This assumption is sufficiently valid if the spatial development of 
the sheet is slow enough for a quasi-two-dimensional situation to 
hold at any downstream location (Reference 2). 

This two-dimensional time-dependent model of Westwater was 
further investigated by many others with the purpose of improving 
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the numerical stability in the calculation. For example, Moore 
(1971) used a point vortex representation to study numerically the 
evolution of an initially plane vortex sheet by using a very 
accurate time integration procedure and two to four times the number 
of point vortices used by Westwater for vortex sheet 
representation. He found that the paths of vortices near the edges 
of the sheet were extremely contorted. This motion is convincingly 
shown to be due to propensity of close line vortices to rotate about 
each other (Reference 3). 

Clement and Maul 1 (1973) used a method similar to that used by 
Westwater to represent the rolling up of a trailing vortex sheet 
(Reference 4). The investigation was oriented towards determining 
the effect of span loading distribution on the roll-up 
characteristics of trailing vortices. 

Hancock (1970) questioned the validity of applying these two- 
dimensional results to finite wings when the lift coefficient is 
large because of the significant effects of wing circulation 
(Reference 5). To improve the two-dimensional model. Butter and 
Hancock (Reference 6) replaced the wing by a single lifting-line 
bound vortex which represented the circulation distribution around 
the wing. The trailing vortex sheet was approximated by a number of 
discrete line vortices, and the downwash and sidewash velocities 
were calculated at midpoint of each vortex element. The induced 
velocity of each vortex element on itself was ignored because the 
self-induced velocity at any point is infinite. A similar approach 
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was also employed by Hackett and Evans for high-lift wings 
(Reference 7). 

In all these investigations, the wing circulation distribution 
was assumed unchanged during the wake roll-up process. In other 
words, the interaction between the wing loading and the wake roll-up 
was not accounted for. The latter can be important for a 
configuration with a canard, or an aft tail, or both, under high 
lift conditions. Numerical methods to determine these interaction 
effects are available and have been applied in the past mostly to 
calculating longitudinal loading on slender wings with edge- 
separated vortex flow (see, for example, References 8 and 9). The 
solution was calculated in an iterative manner by satisfying the 
flow tangency condition on the wing surface and the force-free 
condition on the free vortex sheets. 

In the present three-dimensional study, in addition to 
longitudinal characteristics, the effect of wake shape in sideslip 
and discrete vortices (such as strake or forebody vortex) on lateral 
characteristics will be investigated. The computational method is 
similar to that of Reference 9 in that a lifting-surface method is 
used to represent the wing, not just a lifting line. 

The wing boundary conditions are formulated by the Quasi-Vortex 
Lattice Method (QVLM) of Lan (Reference 10). 

The trailing vortices in the wake region are represented by 
discrete free vortex segments which are aligned with the local 
velocity vector at their midpoint to satisfy the force-free 
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condition as described in detail in Section 2.5.2. The flow 
tangency condition is satisfied on the wing including the trailing 
edge. Due to the nonlinear nature of interaction between the wake 
shape and wing loading, the problem is solved in an iterative 
manner. 

Chapter 2 presents the theoretical method. In Chapter 3, 
numerical results are presented and discussed. Conclusions are made 


in Chapter 4. 



2. THEORETICAL METHOD 


2.1 PROBLEM DESCRIPTION 

For a wing with a positive lift, the sectional loading 
decreases to zero toward the wing tips, producing a circulation 
variation around the wing. This circulation variation causes an 
inward deflection of the streamlines above the wing and an outward 
deflection below the wing. Then the streamlines that converge 
behind the wing are in different directions and form a surface of 
vorticity. This surface of vorticity tends’ to roll up farther 
downstream, forming two distinct vortices of opposite direction. 

The wing can be represented by a horseshoe vortex distribution which 
consists of bound vortices and trailing vortices to create pressure 
differences used to model wing aerodynamic characteristics. The 
wake region is modeled by force-free vortex sheets (continuations of 
wing trailing vortices), across which there is no pressure 
difference. In the present method, the code based on the QVLM 
(Reference 11) that incorporates the wing, fuselage, canard, and 
tail surfaces is modified to determine the wake shape in 
longitudinal and lateral directions and its effects on tail 
planes. It can also determine the effect of canard wake on wing and 
tail planes. 

The following boundary conditions are imposed on the flow 
model : 

a) The flow must be tangent to the wing camber surface. 
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b) The trailing-edge Kutta condition is to be satisfied. 

c) The vortex segments in the wake region are force free. 

In general this nonlinear problem has three sets of unknowns: 

a) strength of the wing horseshoe vortices, 

b) strength of the free vortices in the wake, 

c) the location of wake vortices. 

For simplicity and saving in computing time, the following 
assumptions are made; 

a) The wake vortex strength is assumed equal to wing trailing 
vortex strength. Decay in vortex strength is not 
considered. 

b) The vorticity organizes itself into two or more centers of 
vorticities behind the airplane. This initial roll-up of 
the trailed vorticities behind an airplane is a relatively 
rapid process and generally complete within a distance 
equal to several spans of the wing. For this reason and 
for practical purposes, the part played by viscosity in 
this initial organization is assumed small and the roll-up 
process may, to the first order of accuracy, be computed 
as an inviscid process (Reference 1). 

2.2 GEOMETRY OF WING AND VORTEX WAKE SYSTEM 

In the Quasi-Vortex Lattice Method (Reference 10), for the 
purpose of satisfying wing boundary condition, the continuous vortex 
distribution over the wing is replaced by a quasi-continuous one, 
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being continuous chordwise but stepwise constant in the spanwise 
direction. Thus, the wing surface can be divided into a number of 
vortex strips with the associated trailing vortices (Figure 1). The 
chordwise location of bound elements is given by the cosine law and 
is illustrated in Figure 1 (Reference 10). 


X k = x *e + "f C 1 " cos f( 2k - l)x/2N] } (1) 

k = 1,2 N 

where 

X^ e : leading-edge X-coordinate 

c : chord length 

N : number of bound elements in a chordwise direction. 

A wing may be divided into a number of spanwise sections. In 
each section, the spanwise locations of trailing elements are given 
by the following relation (Reference 10): 



[1 - c° s ( 2 ^ 2 m - )*1 


( 2 ) 


j 


1 , 2 ,... 




where 


M : number of legs of trailing vorticity, which is one 

larger than the number of spanwise strips of bound 
elements 

Y s ^: Y coordinate of the inboard edge of a spanwise 

section 

Y g 2: Y coordinate of the outboard edge of a spanwise 

section 
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The location of control points where the flow tangency 
condition is satisfied is also based on the cosine law, as 
illustrated in Figure 1. The control point locations are as 
follows: 


X 


^ + / ° ' “ s I 1 ’ 


i = 1,. .. ,N 


T sl + H — - - t 1 - C0S( F L) ) 


j = 1,2,. ..(M - 1) 


where 


( 3 ) 


(4) 


X : leading-edge X coordinate at Y_ 

He . c 

1 


j 


Cj : chord length at Y c 


j 


The wake vortex system is made up of wing trailing vortices. 
Each trailed vortex from the trailing edge of the wing is divided 

into a number of segments up to two span lengths behind the wing 

j. 1 

(Figure 2) to form the initial location of the wake vortex system. 
Each segment is assumed to have a length equal to 


AX = .25c 


R 


(5) 


These segments have the following characteristics: 

a) All segments have equal lengths. 

b) The length of each segment is preserved throughout the 
iteration process. 

c) In the final converged position these segments are aligned 
in the direction of the local velocity vector. 
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The initial location of wake vortex segments are 


X = X + iAX 
i te 



where 

L: number of segments in each element. 


( 6 ) 

(7) 

( 8 ) 


2.3 INDUCED VELOCITY DUE TO WING 


As described in Section 2.2, the wing surface is covered with a 
number of vortex strips with the associated trailing vortices. In 
any strip, consider a vortex element of ydX with a direction %. 

The induced velocity due to all bound elements in the i 1 -* 1 strip at 
& = (X, Y, Z) is given by (see Appendix A) 


.2 X te 


V S) - h J u x<x,,( ]yfi Tji l m - ifn* ■ * ,)dx ’ 

and due to the associated trailing vortices by (Reference 10) , 

t 

.2 X„ * ($ - $)xdt 

5 i <ft - h I •)«' 

Y Yl I? 


(9) 


( 10 ) 


£e 
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where symbols a, t>, etc., are position vectors of a vortex element 
relative to a field point, §. = (X, Y, Z) , and are defined in the 
list of Symbols and Appendix A. 
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If the integration variable X' in Equations (9) and (10) is 


replaced by 0, where X' = X + (1 - cose) , then Equations (9) 

£e 2 

and (10) can be reduced to 


$ (£) 

1 


/ ^.(eJYOJsinede 

Sir £ 1 


(ID 


and 


2 77 

($) =-^1^- / $ (0) Y (e)sin0d0 
i 2 8ir 0 2 


(12) 


where 


5 .( 6 ) - - 4-1 ■ *' 

1 | a' x 1' | 2 |S’j 


(13) 



(14) 

(15) 


y(e) = y(x') ; ; 

Let and be $2 for the left and right trailing vortices 
of a bound element, respectively. The total induced velocity due to 
the i tl1 strip of vortex distribution is then given by 


v*> - ^ 


IT 

/ £- ( 0)y( 0)sin0d0 + 

0 


g 2 c(Y) / 

8r l 


S' 2 (e)Y(e)sin9d0 + 


10 



(16) 


— |^ Y) f 5" 2 (9)Y(0)sin0d0 


where the first term is due to bound elements, the second due to the 
left leg of trailing vortices, and the third due to the right leg of 
trailing vortices. The above integrals are reduced to finite sums 
through the midpoint trapezoidal rule (Reference 10) , 


$ ± (&) 


£ 2 c( Y) 
8N 


N 

I (5, + 5' 

k=l k 


+ S-.VinO* 
k 


(17) 


Note that the integration for $ 2 (Equation 14) can be directly 
performed. However, it can also be derived by using Equatiori 
(13). This is indicated in Equation (A. 4) of Appendix A. 


2.4 INDUCED VELOCITY DUE TO WAKE VORTEX SYSTEM 

The wake vortex system, as described in Section 2.2, consists 
of M vortex elements, each of which is divided into a number of 
smaller vortex segments. Assume that the n c ^ element has L small 
segments. The induced velocity at a point (X, Y, Z) , due to the m 1 - 11 
segment of the n 1 -* 1 element is given by (see Appendix A) 




a X t g' 

'|S» X t ’| 2 \t'\ 


V 


where 

.AAA 

R = Xi + Yj + Zk 


(18) 


(19a) 


a - (X - X)i + (Y - Y)j + (Z - Z)k 
m m m 


(19b) 
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(19c) 


s ’ <Vi - X)I + <Vi - + <Vi - z)i 


1 = (X - X )i + (Y - Y )j + (Z - z )k (19d) 

mfl m m+1 m m+1 m 


a' - (X - X)i + g(Y - Y)j + g(Z - Z)k 
m m m 


(19e) 


*' ' ' x>i + s( Vi - Y) j + e< Vi - z)k (I9f> 


*’ ' <Vl - V 1 + 8( Vl - + 8< Vl - V k <19s> 


J.1. 

r : vortex strength of the n tn element, 

n 

The subscripts m and m+1 correspond to the end points of 
the ra 1 -* 1 vortex segment. 

Now, the Induced velocity due to the n th element can be written 
as 


V*> - \ (S) 

m=l m 

Therefore, the induced velocity due to all elements is 

tr(£) = l $ (ft) 

n=l 


( 20 ) 


( 21 ) 


2.5 BOUNDARY CONDITIONS 


As mentioned before, the boundary conditions that are imposed 
on the flow are 

a) The flow must be tangential to the wing camber surface. 

b) The trailing-edge Kutta condition is to be satisfied. 
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c) The vortex elements in the wake behind the trailing edge 
must be force free. ! 


2.5.1 Formulation of Wing Boundary Condition 


The bound elements and the corresponding control points of the 
wing surface are numbered from the leading edge to the trailing edge 
and from root to tip. Thus there are K = N(M - 1) bound elements 
and corresponding control points. The flow tangency condition can 
be written as (References 12 and 13) 

dz c 

[Df j 1 (Yj } = {cos* -r^r- cosa - cos*sina} (22) 


[D* i ] -Cy* j > = (08in* + p (^2 sln< i ) + "kT? cos< f>) 


b/2 


where 


- r 


sin * ( b72 )} 


(23) 


D i.i 


dZ 

c 

dX 


ij 


is the induced normal velocity under symmetrical 
flight conditions at the i*-* 1 control point of 
wing due to a unit horseshoe vortex density at j 
is the vortex density of the j 1 "* 1 bound element 

is the camber slope of the i*^ control point 
is the dihedral angle 
is the induced normal velocity under 
antisymmetrical conditions at the i fc ^ control 
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y' 


j 


3 



00 


point of wing due to a unit horseshoe vortex 
density at j 

is the vortex density of the j 1 -* 1 bound element 
under antisymmetrical conditions 
is the sideslip angle 



OO 


The right hand side of the first equation is the boundary condition 
for symmetrical loading at a given angle of attack. 


2.5.2 Formulation of Vortex Force-Free Condition 

Control points in the wake are chosen at midpoints of vortex 
strips (Figure 2). To satisfy the force-free condition, the vortex 
segment should be aligned in the direction of local velocity vector 
calculated at the control point. For this condition to be 
satisfied, a linear interpolation is used to find the induced 
velocity component at the midpoint of each vortex segment. 

Consider the n c ^ and (n + l)*-* 1 vortex elements and their m 1 -* 1 
vortex segments. The end point coordinates of the m cn vortex 
segments are given by (X n>m , Y n>m , Z n>m ), (X^j, VnrH > Z n,mfl> 

and (Xn+x^m* Y n+1 ,m> z n+l,m^’ ^^n+1, m+l> Y n+1, nrt-1 * Z n+1, m+1^* 

The control point locations are then 
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. (X + X ..) (X +X. 1 _. 1 ) 

v 1 f n,m n,m+l n+l,m n+l,m+l , 

CP = 2 L 2 2 J 


CP 2 


1 r ^n,m + ^n,m+P , ^ Y n+1 ,m + ^n+l,m+l\ 

o" l O O j 


,, 1 r ^ Z n,m + Z n,nri-1^ , ^ Z n+l,m + Z n+l,m+l\ 

Z CP " 2 1 2 + 2 J 


Assume the velocity at the control point of a segment at a 


given iterative step is given by 


$ = (ui + vj + wk) 


where 


u = U cos a + u, 
» 1 


v = v^ , for symmetrical flight condition 


v => v, + U sing , in sideslip 
I 00 


w = U sina + w, 
00 1 


u^ is the induced backwash velocity ' 

Vj^ is the induced sidewash velocity 

w^ is the induced downwash velocity 

U is the free-stream velocity. 

00 

Then, the new location of the (m + l) 1 "* 1 end point will be 


X . - - X + ^ AS 
n,m+l n,m V 


Y , - Y + ^ AS 
n,nri-l n,m V 
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"n.m+l 


Z + £ AS 
n,m V 


where 


(26c) 


, 2 . 2 , 2 
V « /u + v + w 


(27) 


AS = /(X 


n,ra+l 


X ) 2 + (Y 
n,m n,m+l 


Y ) 2 + (Z - Z )' 

n,m n,m+l n,m 


(28) 


Before Equations (26) are used, it should be noted that the 
length of each segment is to be preserved. For numerical stability. 
Equations (26) will be modified as follows by including a relaxation 
parameter. 

Consider the same m tn segment. If this segment moves "a" 
percent (of the initial location) only according to the velocity 
computed at its control point, then Equations (26) can be modified 
to be 


4Y ._ 4S + (1 -a)(Y niWfl 

AZ - Si AS + (1 - a)(Z n ni+1 

o o 2 
AX = /AS - AY - AZ 


" Y ) 
n,m 

- Z ) 
n,m 


where "a" represents a relaxation parameter. 

It follows that 

X . , ** X + AX 
n,m+l n,m 

Y = Y + AY 

n,mrfl n,m 


(29) 


(30a) 

(30b) 
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+ AZ 


(30c) 


Z 


n.uri-1 


Z 

n.m 


2.6 SOLUTION PROCEDURES 


The basic sets of unknowns of the problem are the bound vortex 
density of the wing, and the location of wake vortex elements. The 
problem is nonlinear because the locations of the vortex system are 
unknown a priori. Therefore, the problem will be solved by an 
iterative process described below. 

i) The vortex lattice for the wing surface and the initial 
location for the wake elements are first prescribed. 

ii) By satisfying the wing boundary conditions, the bound 
vortex density of the wing is obtained. 

iii) The free elements of the wake vortex system are adjusted 
to be in the local velocity direction as described in 
Section 2.5.2. 

iv) Steps (ii) through (iii) are repeated until a converged 
solution is obtained. 


After the longitudinal wake shape is obtained, the flow is then 
perturbed by a sideslip angle if lateral characteristics are needed. 

The following procedure is used to find the new asymmetric wake 
shape and aerodynamic forces and moments. 


a) Initially, calculate y (vortex density in sideslip) due 

3 

to gsin<}> using the symmetrical wake shape. 


b) Sum the antisymmetrical and symmetrical vortex densities 


to find a new wake shape under the 0 perturbation through 


the iterative process. 
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c) Calculate the new downwash on the wing with the new 
asymmetrical wake shape. 

d) Subtract the downwash calculated on the wing, with the sum 
of antisymmetrical and symmetrical vortex densities from 
the downwash of symmetrical vortex density to get the 
influence of the new wake. 

e) Add step (d) to 8sin<}> to find new boundary conditions. 

f ) With the new boundary conditions , calculate new vortex 
density. 

g) Calculate all aerodynamic forces and moments. 

The initial locations of the free vortex elements are assumed 
to be in the wing plane. In the iteration process, the force-free 
condition is satisfied on the free elements from the root to tip in 
a row and row by row in the downstream direction. A similar 
approach with different control point locations was used by Butter 
and Hancock (Reference 6) with success. 


2.7 DISCRETE FOREBODY AND STRAKE VORTICES 

Experiments show that for slender bodies and lifting surfaces 
the flow always separates from the surface and rolls up into spiral 
vortex sheets. The effect of this separated-f low or vortex-flow 
phenomenon is quite important at subsonic speed because of its large 
contribution to the total aerodynamic characteristics. 

The leading-edge and side-edge vortex systems induce additional 
velocities at the wing upper surface and produce an additional lift 
force which depends nonlinearly on the angle of attack. 
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In addition, when the leading-edge vortex system passes over an 
area downstream of the leading-edge tip, additional vortex lift is 
developed. This component of lift has been called the augmented 
vortex lift. Also for a straked wing, the strake vortex not only 
produces the augmented vortex lift but also induces downwash inboard 
and upwash outboard of the strake-wing juncture. To account for 
this effect, a discrete vortex starting from the juncture to 
downstream infinity is introduced. The strength of this vortex is 
calculated by equating the vortex lift to the Kutta-Joukowski force 
(Reference 11): 

~ pU OT 2 c s cdY = prw £e d* (31) 

where 

Cg is the sectional leading-edge suction coefficient 
T is the equivalent circulation per unit length of 

the leading edge 

w^^ is the upwash at the leading edge 

d£ is the length measured along the leading edge. 

The total circulation, r , of vortex strength is given as 
follows (Reference 11) 




7 / 


C S C 


, w/U | „ 
£e oo 1 ze 


dY 

d£ 


dZ 


(32) 


where the integration is performed along the leading edge and w is 
evaluated at the leading edge. 
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The average vortex strength, r , per unit length of the leading 
edge Is therefore 



where s^ e Is the length of leading edge on which the vortex is 
produced. 

Note that 


(33) 


if 


b c s c 


2 J w/U I 
a “>' Ze 


1 1 ^ c< ? { ' 
dY = j [-kb - a) / 


0 w/U J*e 


sin<j>d(j)] 


■*<>-•> 5.1 


CgC 

. , w/U I 8in<l, i 
i= >1 <»' Ze 


(34) 


b) Since both c c and w„ depend on , the circulation 
s Ze t 

is determined iteratively. 

A similar approach to that in Section 2.5.2 is used to find the 
position of the discrete vortex element from the wing junction to 
distance of two spans behind the wing trailing edge. Also the 

t 

position of discrete vortex element for small sideslip perturbation 
can be calculated. 

The main objective of finding the location of discrete vortex 
elements is to determine the location of this vortex relative to the 
vertical tail, and to determine its effect on directional 
stability. Additional discussion will be given in the next chapter. 
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3. NUMERICAL RESULTS AND DISCUSSIONS 


The present numerical method, which is based on discretization 
of the vortex sheet into vortex segments, is used to verify the 
roll-up process of the trailing vortices. Also for any surface with 
a kinked leading edge, the trailing vortices from the kinked 
location is calculated and demonstrated clearly. 

The following points should be considered in applications of 
the method: 

a) The accuracy of the method for determining longitudinal 
and lateral-directional derivatives and wake shape in' 
sideslip depends somewhat on the length of each free 
vortex segment and the number of trailing vortices. If 
the length of free vortex segment decreases , the method 
tends to be more accurate; but the computing time will be 
increased. 

b) The complete trailing vortex system is calculated by the 
step-by-step process described in Chapter 2. Generally 
the roll-up is completed between the fifth and tenth 
iterations. 

The longitudinal wake shape by the present method is compared 
with the numerical solution of Bloom and Jen, where they applied the 
artificial viscosity method of Kuwahara and Takami to a rectangular 
wing with aspect ratio of eight and the lift coefficient equal to 
one (C L = 1; Reference 14). The results of the present method were 
calculated by using 15 spanwise vortex strips and each trailing 
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element divided into twenty segments. The length of each free 
vortex segment was 25 percent of the root chord. It took 
approximately 10 iterations for the complete roll-up process. The 
predicted wake shape is plotted in Figure 3, which shows a good 
agreement between the two methods. 

The longitudinal roll-up process of trailing vortices of a 
twin- jet fighter airplane, at two different angles of attack, from 
the wing trailing edge up to two spans behind the trailing edge, is 
calculated. The fuselage is represented by a body of revolution, 
and nacelles are not modeled. The results indicate that at low 
angles of attack, the roll-up process is slow because of small 
loading as well as weak interaction among trailing vortices (Figure 
4). At high angles of attack, the loading is increased and the 
roll-up process is much faster as is shown in Figure 5. For a = 20 
degrees, the roll-up is completed within two spans behind the 
trailing edge. The rapid change In wake shape at Y/(b/2) « 0.35 in 
Figure 5(g) is due to change of dihedral angle in the tip section 
and discontinuity in the wing leading edge. This point will move 
inboard as the roll-up process is completed. 

The aerodynamic characteristics of the twin-jet fighter 
airplane configuration are compared with the experimental data of 
Reference 15 and QVLM (Reference 11; Inviscid, flat wake assumption, 
and with vortex breakdown) in Figure 6. The results are in good 
agreement up to an angle of attack of about 15 degrees. The 
disagreement with experimental data above 15 degrees is perhaps due 
to viscous effects. 
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The predicted longitudinal aerodynamic characteristics of tail 
plane by the present method is compared with the QVLM results in 
Figure 7. The results show some disagreement at high angles of 
attack. This is because, with roll-up, the wake tends to move away 
from the low tail at high angles of attack so that the induced 
downwash is reduced. 

The effect of wake in sideslip for the twin-jet fighter 
airplane is calculated and the results are shown in Figure 8. It is 
seen that the wake is not symmetrical any more and moves in the 
direction of sideslip. 

The predicted lateral-directional characteristics of the twin- 
jet airplane in sideslip are compared with the QVLM results 
(inviscid, flat wake, and with vortex breakdown) and experimental 
data of Reference 16 in Figure 9. The present method shows a small 
deviation from the QVLM results , which is due to the effect of wake 
roll-up on the vertical tail. 

The longitudinal and lateral wake shape of a straked-wing 
configuration in sideslip is determined and the results are given in 
Figures 10-11. It appears that the strake vortex tends to move 
outboard and stay close to the tip vortex as it moves downstream. 

The predicted longitudinal and lateral-directional aerodynamic 
characteristics of a straked-wing configuration in sideslip are 
compared with the QVLM results (inviscid, flat wake, and without 
vortex breakdown) and experimental data of Reference 17 in Figures 
12-14. Significant changes in longitudinal aerodynamic 
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characteristics due to tail plane and lateral-directional stability 
derivatives are observed as follows: 

a) The longitudinal vortex wake of the straked-wing 
configuration is rolled up within a distance of less than 
one span behind the trailing edge. The presence of the 
strake vortex introduces some effect on the downwash field 
at the horizontal tail. For this reason, some changes in 
aerodynamic characteristics of the tail plane calculated 
by the present method relative to the QVLM results 
(inviscid, flat wake, and without vortex breakdown) are 
observed (Figure 12). The and results from QVLM 
and the present method for the complete configuration, 
show good agreement with experimental data (without 
horizontal tail) of Reference 17 at low angles of 
attack. At high angles of attack, because of inviscid 
assumption, the above methods failed to predict the 
aerodynamic characteristics accurately (Figure 13(a)). 

Due to change of downwash field of the tall plane by the 
roll-up wake, the predicted pitching moment coefficient 
(C m ) by the present method is slightly better than the 
QVLM prediction (inviscid, flat wake assumption, and 
without vortex breakdown) as shown in Figure 13(b). 

b) The predicted C m is more negative than the data show. The 
calculated results indicate that this is mostly due to the 
contribution of augmented vortex lift. Since the latter 
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is assumed to be uniformly distributed along the wing- 
strake junction, further research on its distribution is 
needed to improve the modeling. 

c) Due to the presence of asymmetrical vortex wake roll-up of 

straked wing in sideslip, the flow pattern around the 

vertical tail changes and introduces a side force which 

causes significant changes in directional stability 

derivative (C ) and also c derivative of the 
n„ Y„ 

0 0 

airplane. Therefore for airplanes with high wing 
loadings, the assumption of flat wake in calculating the 
lateral-directional derivatives should be reexamined 
(Figure 14). 

Figures 15 and 16 show the position of discrete strake vortices 
for the straked-wing configuration at the vertical tail location. 

No numerical or experimental data are available for comparison. 
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4. CONCLUSIONS 


A theoretical method has been developed for predicting the wake 
shapes in symmetrical flight conditions and with sideslip, including 
their effect on tail planes. The present method has been shown to 
work satisfactorily for most configurations without encountering any 
numerical instability in the solution. The results indicate that 
for a wing with high lift, the assumption of flat wake is 
inaccurate; and the effect of vortex wake roll-up on longitudinal 
and lateral-directional aerodynamic characteristics of airplanes 
needs to be examined. Note that in the present method the effect of 
fuselage vortex wake roll-up was not accounted for. 
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APPENDIX A: EVALUATION OF INDUCED VELOCITY DUE TO A HORSESHOE 

VORTEX AND LINE VORTEX SEGMENT 


The velocity field induced by a line vortex element of strength 
T is given by (Reference 10) 


$(&) = 


e 2 r 

4ir ' 


(£, - tl) x d£ 

/ — 5 

1 E S 


where 

0 = /I - M 2 

■ A A A 

R = Xi + Yj + Zk 


(A. 1) 


■ A A A 

R, = X'i + Y'j + Z’k 
1 J 

R 2 o ■ (X - X') 2 + 0 2 (Y - Y') 2 + 0 2 (Z - Z') 2 

3 

where (X', Y' , Z') is the point on the vortex element. 

The horseshoe vortex consists of a bound element and two 
trailing vortices (Figure 17). The induced velocity due to the 
bound element with end points (X^, Y^ , ) and (X£, Y^ , Z^), as 

indicated in Figure 17, is obtained by substituting 

& = a + rl , with t representing a fraction of the vortex 
element, into Equation (A.l) as follows: 
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2(2A + 6) 

(B 2 - 4AC)(A + B + C) 


172 


> , B 2 - 4AC ± 0 


(A. 2) 


where 

A = |t'| 2 . B = 2a' • P , C = |a* | 2 • 
Further it can be shown that (Reference 10) 


B 2 - 4AC = -4|a’ x V \ 2 
2k + B = 2$’ • £' 

A + B + C = |S' | 2 

where 

a = (X t - X)i + (Y : - Y)j + (Z l - Z)k 

t = (X 2 - X) i + (Y 2 - Y) j + (Z 2 - Z)k 

1 = (X 2 - Xj)! + (Y 2 - Yj)j + (Z 2 - Z L )k 

a* = (X 1 - X)i + B(Y X - Y) j + g(Z 1 - Z)k 

; I . 

b' = (X 2 - X)i + B(Y 2 - Y)j + 0(Z 2 - Z)k 

V = <X 2 - X 1 )i + B(Y 2 - Y)j + B(Z 2 - Z^k 

On rearranging. Equation (A. 2) becomes 


5 < s > ■ 17 


a x 1 


xt'T |£' 


-} • V 


(A. 3) 
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For the right trailing vortex element. Equation (A. 3) becomes 


(if the trailing legs go to infinity) 


2_ x i 


_ a' • i 2 

V ft) - ,7; T72 {1 T^r } ~ 4* G 2 


(A. 4) 


|a’ 2 x i | 


a\ 


where 


a 2 = (x 2 - X)i + (Y 2 - Y)j + (Z 2 - Z)k 


a’ 2 « - X)i + g(Y 2 - Y)j + g(Z 2 - Z)K 


For the left trailing vortex of a bound element. Equation (A. 4) is 
still applicable except that a negative sign is needed and the 
subscript ”2" in a and a' is replaced with "1”. 

Note that for the induced velocity of free vortex segments , 
Equation (A. 3) is also applicable. 
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APPENDIX B: GEOMETRY DATA FOR EXAMPLE CONFIGURATIONS 


Three different configurations used in this study are as 
follows: 

1. Rectangular wing (Reference 14) 

2. Twin-jet fighter airplane (Reference 16) 

3. Straked wing (Reference 17) 

The geometry of the above configurations is presented as 
follows : 

a) For the rectangular wing, see Figure 18. 

b) For the twin-jet fighter airplane, see Table I. 

c) For the straked wing configuration, refer to Figure 19 
and Table II. 
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table I 


Twin-Jet Fighter Airplane Geometry (Taken from Ref.' 16) 

1 DIMENSIONAL CHARACTERISTICS OF AIRPLANE 


Overall length . . . 

Wing: 

Span . . 

Area (Including leading-edge extension) . . . 

Root chord 

Tip chord 

Mean aerodynamic chord, c 

Leading edge of c rearward of leading edge 

of root chord . ■ 

Aspect ratio 

Taper ratio 

Sweepback of 25-per.cent chord 

Dihedral (Inboard 69.5 percent b/2) . . . . . 
Dihedral (outboard 69.5 percent b/2) . . . . 

Incidence . 

Airfoil section: 

Root 

Tip . . 

Aileron: 

Area (one side) rearward of hinge line ... 
Span (one aileron) (from 44.5 to 

67.0 percent b/2) 

Inboard end chord (base line 

103.24 in. (262.23 cm)) 

Outboard end chord (base line 

155.44 In. (394.82 cm)) 

Spoilers: 

Area (one side) 

Span (from 45.3 to 67.0 percent b/2) . . , 

Inboard end chord 

Outboard end chord 

Horizontal tail: 

Area (in chord plane) , 

Movable area 

Span . . . , 

Mean aerodynamic chord of horizontal tail 
Aspect ratio 


57.59 ft (17.55 m) 

38.41 ft (11.71 m) 

538.34 ft 2 (50.01 m2) 

282.00 in. (716.28 cm) 

47.00 in. (119.38 cm) 

192.50 in. (488.95 cm) 

........ 110.76 in. (281.33 cm) 

2.82 

. . . 0.167 

45.00° 

0° 

12.00° 

1.00° 

NACA 0006.4-64 (modified) 

NACA 0003.0-64 (modified) 

13.08 ft 2 (1.22 m2) 

. 4.35 ft (1.33 m) (22.5 percent b/2) 

37.81 in. (96.04 cm) (21.3 percent Z) 

34.38 in. (87.33 cm) (27.6 percent Z) 

: 5.44 ft 2 (0.50 m 2 ) 

4.19 ft (1.28 m) 

1.39 ft (0.42 m) 

1.04 ft (0.32 m) 

94.70 ft 2 (8.80 m 2 ) 

77.40 ft 2 (7.19 m 2 ) 

17.705 ft (5.40 m) 

6.143 ft (1.87 m) 

3.30 
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TABLE I 


Twin-Jet Fighter Airplane Geometry (Taken from Ref. 16 ) 
DIMENSIONAL CHARACTERISTICS OF AIRPLANE - Concluded 


Taper ratio 

Sweepback of 25-percent chord 

Dihedral 

Root chord (at airplane center line) 

Tip chord (theoretical) 

Airfoil section: 

Root (airplane center line) 

Tip (theoretical) 

Hinge-line location, percent c^ ........ . 

Vertical tail: 

Area 

Span 

Taper ratio 

Root chord 

Tip chord 

Sweepback of 25-percent chord 

Airfoil section: 

Root 

Tip 

Rudder: 

Area (rearward of hinge line) 

Hinge-line location, percent of water-line chords 


0.20 

35.50° 

-23.00° 

. . 107.00 in. (271.78 cm) 
. . . 21.40 in. (54.36 cm) 

NACA 0003.7-64 (modified) 
NACA 0003.0-64 (modified) 
41.00 

. . . 67.50 ft2 (6.27 m2) 

6.38 ft (1.94 m) 

0.227 

. . 207.15 in. (526.16 cm) 
. . 47.10 in. (119.63 cm) 
58.30° 

NACA 0004.0-64 (modified) 
NACA 0002.5-64 (modified) 

. . . .11.07 ft2 (1.03 m2) 
80.00 
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TABLE II 


Strake-Wing Geometry (Taken from Ref. 17) 



Y 


7v 

AD iq ■ — - — /L+. 

X 

y 

cm 

in. 

cm. 

in . 

0.000 

0.000 

0.000 

0.000 

0.455 

0.179 

0.229 

0.090 

1.115 

0.439 

0.457 

0.180 

1.935 

O.762 

0.68 6 

0.270 

2.893 

1.139 

0.914 

0.360 

3.980 

1.567 

1.143 

0.450 

5.189 

2.043 

1.374 

0.541 

7.968 

3.137 

1.831 

0.721 

11.247 

4.428 

2.289 

0.901 

15.085 

5.939 

2.746 


19.558 

7.700 

3.203 

I.261 

24.491 

9.642 

3.660 

1.441 

30.170 

11 . 878 

4.120 

1.622 

38.892 

15.312 

4.577 

1.802 
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Fig. 2 
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Fig. 4(a) The Wake Shape of a Wing-Body-Tail 

i 

Configuration at Alpha = 0.5 deg. 
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Fig. 4(b) The Wake Shape of a Wing-Body-Tail Configuration 

at Alpha = 0.5 deg. 
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Fig. 4(c) The Wake Shape of a Wing-Body-Tail Configuration 
at Alpha = 0.5 deg. 


41 



2/(b/2) Z/(b/2) 


Fig. 4(d) The Wake Shape of a Uing-Body-Tail Configuration 

at Alpha =0.5 deg. 
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Fig. 4(e) The Wake Shape of a Wing-Body-Tail Configuration 

at Alpha = 0.5 deg. 
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Fig. 5(a) The Wake Shape of a Wing-Body-Tail Configuration 

at Alpha = 20.0 deg. 
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Fig. 5(b) The Wake Shape of a Wing-Body-Tail Configuration 

at Alpha = 20.0 deg. 
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Z/(b/2.) 


□—Woke Shcpe at X=2.80*(b/2.) (m) 
O-Woke Shape at X=3.11«(b/2.) (m) 



Y/(b/2.) 

Fig. 5(c) The Wake Shape of a Wing-Body-Tail Configuration 

at Alpha = 20.0 deg. 
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Z/(b/2.) 


□— Woko Shape at X=3.41 •(b/2.) (m) 
O-Wako Shape at X-3.71«(b/2.) (m) 



at Alpha = 20.0 deg. 
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□ -Wake Shape at X=4.02*(b/2.) (m) 
Q-Wake Shape at X=4.32*(b/2.) (m) 
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Fig. 5(e) The Wake Shape of a Wing-Body-Tail Configuration 

at Alpha = 20.0 deg. 
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Fig. 5(f) The Wake Shape of a Wing-Body-Tail Configuration 

at Alpha = 20.0 deg. 
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Fig. 5(g) The Wake Shape of a Wing-Body-Tail Configuration 

at Alpha = 20.0 dee. 
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Lift Coefficient, CL Drag Coefficient, CD 




Fig. 6(a) Static Longitudinal Characteristics of a 
Wing-Body-Tail Configuration 
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Pitching Moment Coefficient, CM 



0.0 5.0 10.0 15.0 20.0 25.0 30.0 

Angle of Attack, Alpha (deg) 

Fig. 6(b) Static Longitudinal Characteristics of a 
Wing-Body-Tail Configuration (Cont) 
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Horizontal Tall lift Coefficient, CL 



Fig. 7(a) Static Longitudinal Characteristics of Tail Plane of a 
Wing-Body-Tail Configuration 
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Horizontal Tati Drag Coefficient, CD 



Fig. 7(b) Static Longitudinal Characteristics of Tail Plane 
of a Wing -Body-Tail Configuration 
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Horizontal Tall Pitching Moment Coefficient, CM Q 33 ~ 



Fig. 7(c) Static Longitudinal Characteristics of Tail Plane 
of a Wing-Body-Tail Configuration 
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Fig. 8(a) The Wake Shape of a Wing_Body_Tail Configuration 
at Alpha = 20 and Beta = 5 deg. 




□ —Right Wing Wake Shape at X=2.3*(b/2) 
©-Left Wing Wake Shape at X=2.3*(b/2) 
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Fig. 8(b) The Wake Shape of a Wing_Body_Tail Configuration 
at Alpha = 20 and Beta = 5 deg. 




Fig. 8(c) The Wake Shape of a Wing_Body_Tail Configuration 
at Alpha = 20 and Beta = 5 deg. 




Fig. 8(d) The Wake Shape of a Wing_Body_Tail Configuration 
at Alpha = 20 and Beta = 5 deg. 




Fig. 8(e) The Wake Shape of a Wing_Body_Tail Configuration 
at Alpha = 20 and Beta = 5 deg. 





Fig. 8(f) The Wake Shape of a Wing_Body_Tail Configuration 
at Alpha = 20 and Beta = 5 deg. 



□—Right Wing Wake Shape at X=3.5*(b/2) 
G-Left Wing Wake Shape at X=3.5*(b/2) 
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Fig. 8(g) The Wake Shape of a Wing_Body_Tail Configuration 
at Alpha = 20 and Beta = 5 deg. 




Fig. 8(h) The Wake Shape of a Wing _ Body Tall Configuration 

at Alpha = 20 and Beta = 5 deg. 




Fig. 8(i) The Wake Shape of a Wing_Body_Tail Configuration 
at Alpha = 20 and Beta = 5 deg. 




Fig* 80) The Wake Shape of a Wing_Body_Tail Configuration 
at Alpha = 20 and Beta = 5 deg. 
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Fig. 8(k) The Wake Shape of a Wing_Body_Tail Configuration 
at Alpha = 20 and Beta = 5 deg. 
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Fig. 8(1) The Wake Shape of a Wing_Body_Tail Configuration 
at Alpha = 20 and Beta = 5 deg. 




Fig. 8(m) The Wake Shape of a Wing_Body_Tail Configuration 
at Alpha = 20 and Beta = 5 deg. 
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Fig. 8(n) The Wake Shape of a Wing_Body_Tail Configuration 
at Alpha = 20 and Beta = 5 deg. 
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Fig. 8(o) The Wake Shape of a Wing_Body_Tail Configuration 
at Alpha = 20 and Beta = 5 deg. 
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Fig. 9(a) Static Lateral— Directional Stability Derivatives 
of a Wing_Body_Tail Configuration 
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Inviscid Theory with Wake Deformation 
Inviscid Theory without Wake Deformation 

Experimental Data (Ref. 16 ) 


Angle of Attack, Alpha (deg) 


Fig. 9(b) Static Lateral— Directional Stability Derivatives 
of a Wing_Body_Tail Configuration 
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Fig. 10(a) Wake Shape of the Straked Wing-Body Configuration 

at Alpha = 29. deg. 
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Fig. 10(b) Wake Shape of the Stroked Wing-Body Configuration 

at Alpha = 29. deg. 
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Fig. 10(c) Wake Shape of the Straked Wing-Body Configuration 

at Alpha = 29. deg. 
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Fig. 10(d) Wake Shape of the Straked Wing-Body Configuration 

at Alpha = 29. deg. 
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Fig. 10(e) Wake Shape of the Straked Wing-Body Configuration 

at Alpha = 29. deg. 
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Fig. 10(f) Wake Shape of the Straked Wing-Body Configuration 

at Alpha = 29. deg. 


79 




0.0 0.5 1.0 1.5 


Y/(b/2) 

Fig. 10(g) Wake Shape of the Straked Wing-Body Configuration 

at Alpha = 29. deg. 
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-Wako Shape at X=6.6*(b/2) 
•Wake Shape at X=6.8*(b/2) 



Fig. 10(h) Wake Shape of the Straked Wing-Body Configuration 

at Alpha = 29, deg. 
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Right Wing Wake Shape at X«=3.3*(b/2) 
Left wing Wake Shape at X-3.3*(b/2) 
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Fig, 11(a) Wake Shape of the Straked Wing -Body Configuration 
at Alpha = 29 and Beta =17 deg. 





Right Wing Wake Shape at X=3.5*(b/2) 
Left wing Wake Shape at X-3.5*(b/2) 
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Fig. 11(b) Wake Shape of the Straked Wing-Body Configuration 
at Alpha = 29 and Beta = 17 deg. 




-Right Wing Wake Shape at X=3.7*(b/2) 
•Left wing Wake Shape at X-3.7*(b/2) 
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Fig. 11(c) Wake Shape of the Straked Wing-Body Configuration 
at Alpha = 29 and Beta = 17 deg. 



Right Wing Wake Shape at X=3.9*(b/2) 
Left wing Wake Shape at X-3.9*(b/2) 
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Fig. 11(d) Wake Shape of the Straked Wing -Body Configuration 


at Alpha = 29 and Beta =17 deg. 




Right Wing Waka Shape at X«4.2*(b/2) 
Laft wing Waka Shapa at X-4.2*(b/2) 
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Right Wing Wake Shape at X=*4.4*(b/2) 
Left wing Wake Shape at X«4.4*(b/2) 
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Fig. 11(f) Wake Shape of the Straked Wing-Body Configuration 
at Alpha = 29 and Beta = 17 deg. 




Right Wing Wake Shape at X«4.7*(b/2) 
Left wing Wake Shape at X-4.7*(b/2} 
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Fig. 11(g) Wake Shape of the Straked Wing-Body Configuration 


at Alpha = 29 and Beta = 17 deg. 



Right Wing Woke Shape at XM.9*(b/2) 
Left wing Wake Shape at X—4.9*(b/2) 
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Fig. 11(h) Wake Shape of the Straked Wing-Body Configuration 
at Alpha = 29 and Beta = 17 deg. 



Right Wing Wake Shape at X*=5.2*(b/2) 
Left wing Wake Shape at X-5.2*(b/2) 


Z/(b/2) 



Y/(b/2) 


Fig. 11 (i) Wake Shape of the Straked Wing -Body Configuration 
at Alpha = 29 and Beta = 17 deg. 




-Right Wing Woke Shop® at X=5.4*(b/2) 
•Left wing Wake Shape at X~5.4*(b/2) 
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Fig. ll(j) Wake Shape of the Straked Wing-Body Configuration 
at Alpha = 29 and Beta = 17 deg. 




Right Wing Woke Shape at X=5.7*(b/2) 
Left wing Wake Shape at X— 5.7*(b/2) 
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Fig. ll(k) Wake Shape of the Stral 


at Alpha = 29 and 





Right Wing Woke Shape at X=5.9*(b/2) 
Left wing Wake Shape at X— 5.9*(b/2) 
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Fig. 11(1) Wake Shape of the Straked Wing-Body Configuration 


at Alpha = 29 and Beta = 17 deg. 



Right Wing Wake Shape at X*=6.1*(b/2) 
Left wing Wake Shape at X-6.1*(b/2) 
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Fig* 11 (ro) Wake Shape of the Straked Wing— Body Configuration 


at Alpha = 29 and Beta = 17 deg. 






Right Wing Woke Shape at X=6.6*(b/2) 
Left wing Wake Shape at X-6.6*(b/2) 
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Fig. 11 (o) Wake Shape of the Strak 


at Alpha = 29 and 






Right Wing Wake Shape at X*=6.9*(b/2) 
Left wing Wake Shape at X«*6.7»(b/2) 
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Fig. 11 (p) Wake Shape of the Straked Wing-Body Configuration 
at Alpha = 29 and Beta = 17 deg. 



5 


-Inviscid Theory with Wake Deformation 
. Inviscid Theory without Wake Deformation 

Baaed on Wing Reference Area 
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Inviscid Theory with Wake Deformation 

-Inviscid Theory without Wake Deformation 

CD t Baaed on Wing Reference Area 



Angle of Attack. Alpha (deg) 


Fig. 12(b) Static Longitudinal Characteristics of the 


Straked Wing-Body Configuration 
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Inviscid Theory with Wake Deformation 

. Inviscid Theory without Wake Deformation 

CM Based on Wing Geometry 
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Fig. 13(a) Static Longitudinal Characteristics of the 

Straked Wing-Body Configuration at Beta = 0.0 deg. 
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Fig. 13(b) Static Longitudinal Characteristics of the 
Straked Wing-Body Configuration (Cont) 


102 


Effective -Dihedral Derivative, CLB Directional-Stability Derivative, CNB Side-Force Derivative, 


- 0 . 02 1 - 


- 0.01 


0.00 


-(rad -1 ) 

- 

Inviscid Theory 

with Wake Deformation 
without Wake Deformation 

i L 

. 

1 — i 1 . 1 

J 1 L. 1 


10 


15 


20 


25 


20 




Fig. 14 Static Lateral-Directional Stability Derivatives 
of the Straked Wing-Body Configuration at Beta =* 17 deg. 
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Q — Stroke Vortex Position 
Vertical Tail 
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15 Effect of Strake Vortex on Vertical Tail for Straked 
Wing-Body Configuration at Alpha = 29 deg. 
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Fig. 16 Effect of Strake Vortex on Vertical Tail for Straked 
Wing-Body Configuration at Alpha = 29 and Beta = 17 deg. 
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• Span (b) =8 
Chord Length (c) = 1 
Lift Coefficient = 1 

Rectangular Wing Geometry and Aerodynamic Configuration 

Fig. 18 
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